Energy metabolism modulates the regulatory impact of activators on gene expression.

Gene expression is a regulated process fueled by ATP consumption. Therefore, regulation must be coupled to constraints imposed by the level of energy metabolism. Here, we explore this relationship both theoretically and experimentally. A stylized mathematical model predicts that activators of gene expression have variable impact depending on metabolic rate. Activators become less essential when metabolic rate is reduced and more essential when metabolic rate is enhanced. We find that, in the Drosophila eye, expression dynamics of the yan gene are less affected by loss of EGFR-mediated activation when metabolism is reduced, and the opposite effect is seen when metabolism is enhanced. The effects are also seen at the level of pattern regularity in the adult eye, where loss of EGFR-mediated activation is mitigated by lower metabolism. We propose that gene activation is tuned by energy metabolism to allow for faithful expression dynamics in the face of variable metabolic conditions.

Energy metabolism modulates the regulatory impact of activators on gene expression.

Gene expression is a regulated process fueled by ATP consumption. Therefore, regulation must be coupled to constraints imposed by the level of energy metabolism. Here, we explore this relationship both theoretically and experimentally. A stylized mathematical model predicts that activators of gene expression have variable impact depending on metabolic rate. Activators become less essential when metabolic rate is reduced and more essential when metabolic rate is enhanced. We find that in the Drosophila eye, expression dynamics of the yan gene are less affected by loss of EGFR-mediated activation when metabolism is reduced, and the opposite effect is seen when metabolism is enhanced. The effects are also seen at the level of pattern regularity in the adult eye, where loss of EGFR-mediated activation is mitigated by lower metabolism. We propose that gene activation is tuned by energy metabolism to allow for faithful expression dynamics in the face of variable metabolic conditions.

Ratiometric sensing of Pnt and Yan transcription factor levels confers ultrasensitivity to photoreceptor fate transitions in Drosophila.

Cell state transitions are often triggered by large changes in the concentrations of transcription factors and therefore large differences in their stoichiometric ratios. Whether cells can elicit transitions using modest changes in the ratios of co-expressed factors is unclear. Here, we investigate how cells in the Drosophila eye resolve state transitions by quantifying the expression dynamics of the ETS transcription factors Pnt and Yan. Eye progenitor cells maintain a relatively constant ratio of Pnt/Yan protein, despite expressing both proteins with pulsatile dynamics. A rapid and sustained twofold increase in the Pnt/Yan ratio accompanies transitions to photoreceptor fates. Genetic perturbations that modestly disrupt the Pnt/Yan ratio produce fate transition defects consistent with the hypothesis that transitions are normally driven by a twofold shift in the ratio. A biophysical model based on cooperative Yan-DNA binding coupled with non-cooperative Pnt-DNA binding illustrates how twofold ratio changes could generate ultrasensitive changes in target gene transcription to drive fate transitions. Thus, coupling cell state transitions to the Pnt/Yan ratio sensitizes the system to modest fold-changes, conferring robustness and ultrasensitivity to the developmental program.

A new approach for extracting information from protein dynamics.

Increased ability to predict protein structures is moving research focus towards understanding protein dynamics. A promising approach is to represent protein dynamics through networks and take advantage of well-developed methods from network science. Most studies build protein dynamics networks from correlation measures, an approach that only works under very specific conditions, instead of the more robust inverse approach. Thus, we apply the inverse approach to the dynamics of protein dihedral angles, a system of internal coordinates, to avoid structural alignment. Using the well-characterized adhesion protein, FimH, we show that our method identifies networks that are physically interpretable, robust, and relevant to the allosteric pathway sites. We further use our approach to detect dynamical differences, despite structural similarity, for Siglec-8 in the immune system, and the SARS-CoV-2 spike protein. Our study demonstrates that using the inverse approach to extract a network from protein dynamics yields important biophysical insights.

A new approach for extracting information from protein dynamics.

Increased ability to predict protein structures is moving research focus towards understanding protein dynamics. A promising approach is to represent protein dynamics through networks and take advantage of well-developed methods from network science. Most studies build protein dynamics networks from correlation measures, an approach that only works under very specific conditions, instead of the more robust inverse approach. Thus, we apply the inverse approach to the dynamics of protein dihedral angles, a system of internal coordinates, to avoid structural alignment. Using the well-characterized adhesion protein, FimH, we show that our method identifies networks that are physically interpretable, robust, and relevant to the allosteric pathway sites. We further use our approach to detect dynamical differences, despite structural similarity, for Siglec-8 in the immune system, and the SARS-CoV-2 spike protein. Our study demonstrates that using the inverse approach to extract a network from protein dynamics yields important biophysical insights.

Fly-QMA: Automated analysis of mosaic imaginal discs in Drosophila.

Mosaic analysis provides a means to probe developmental processes in situ by generating loss-of-function mutants within otherwise wildtype tissues. Combining these techniques with quantitative microscopy enables researchers to rigorously compare RNA or protein expression across the resultant clones. However, visual inspection of mosaic tissues remains common in the literature because quantification demands considerable labor and computational expertise. Practitioners must segment cell membranes or cell nuclei from a tissue and annotate the clones before their data are suitable for analysis. Here, we introduce Fly-QMA, a computational framework that automates each of these tasks for confocal microscopy images of Drosophila imaginal discs. The framework includes an unsupervised annotation algorithm that incorporates spatial context to inform the genetic identity of each cell. We use a combination of real and synthetic validation data to survey the performance of the annotation algorithm across a broad range of conditions. By contributing our framework to the open-source software ecosystem, we aim to contribute to the current move toward automated quantitative analysis among developmental biologists.

Large-scale analysis of micro-level citation patterns reveals nuanced selection criteria.

The analysis of citations to scientific publications has become a tool that is used in the evaluation of a researcher's work; especially in the face of an ever-increasing production volume1-6. Despite the acknowledged shortcomings of citation analysis and the ongoing debate on the meaning of citations7,8, citations are still primarily viewed as endorsements and as indicators of the influence of the cited reference, regardless of the context of the citation. However, only recently has attention9,10 been given to the connection between contextual information and the success of citing and cited papers, primarily because of the lack of extensive databases that cover both types of metadata. Here we address this issue by studying the usage of citations throughout the full text of 156,558 articles published by the Public Library of Science (PLoS), and by tracing their bibliometric history from among 60 million records obtained from the Web of Science. We find universal patterns of variation in the usage of citations across paper sections11. Notably, we find differences in microlevel citation patterns that were dependent on the ultimate impact of the citing paper itself; publications from high-impact groups tend to cite younger references, as well as more very young and better-cited references. Our study provides a quantitative approach to addressing the long-standing issue that not all citations count the same.

Single-Cell Transcriptomic Analysis of Human Lung Provides Insights into the Pathobiology of Pulmonary Fibrosis.

Rationale: The contributions of diverse cell populations in the human lung to pulmonary fibrosis pathogenesis are poorly understood. Single-cell RNA sequencing can reveal changes within individual cell populations during pulmonary fibrosis that are important for disease pathogenesis. Objectives: To determine whether single-cell RNA sequencing can reveal disease-related heterogeneity within alveolar macrophages, epithelial cells, or other cell types in lung tissue from subjects with pulmonary fibrosis compared with control subjects. Methods: We performed single-cell RNA sequencing on lung tissue obtained from eight transplant donors and eight recipients with pulmonary fibrosis and on one bronchoscopic cryobiospy sample from a patient with idiopathic pulmonary fibrosis. We validated these data using in situ RNA hybridization, immunohistochemistry, and bulk RNA-sequencing on flow-sorted cells from 22 additional subjects. Measurements and Main Results: We identified a distinct, novel population of profibrotic alveolar macrophages exclusively in patients with fibrosis. Within epithelial cells, the expression of genes involved in Wnt secretion and response was restricted to nonoverlapping cells. We identified rare cell populations including airway stem cells and senescent cells emerging during pulmonary fibrosis. We developed a web-based tool to explore these data. Conclusions: We generated a single-cell atlas of pulmonary fibrosis. Using this atlas, we demonstrated heterogeneity within alveolar macrophages and epithelial cells from subjects with pulmonary fibrosis. These results support the feasibility of discovery-based approaches using next-generation sequencing technologies to identify signaling pathways for targeting in the development of personalized therapies for patients with pulmonary fibrosis.

Within-Gene Shine-Dalgarno Sequences Are Not Selected for Function.

The Shine-Dalgarno (SD) sequence motif facilitates translation initiation and is frequently found upstream of bacterial start codons. However, thousands of instances of this motif occur throughout the middle of protein coding genes in a typical bacterial genome. Here, we use comparative evolutionary analysis to test whether SD sequences located within genes are functionally constrained. We measure the conservation of SD sequences across Enterobacteriales, and find that they are significantly less conserved than expected. Further, the strongest SD sequences are the least conserved whereas we find evidence of conservation for the weakest possible SD sequences given amino acid constraints. Our findings indicate that most SD sequences within genes are likely to be deleterious and removed via selection. To illustrate the origin of these deleterious costs, we show that ATG start codons are significantly depleted downstream of SD sequences within genes, highlighting the constraint that these sequences impose on the surrounding nucleotides to minimize the potential for erroneous translation initiation.

A novel framework for evaluating the performance of codon usage bias metrics.

The unequal utilization of synonymous codons affects numerous cellular processes including translation rates, protein folding and mRNA degradation. In order to understand the biological impact of variable codon usage bias (CUB) between genes and genomes, it is crucial to be able to accurately measure CUB for a given sequence. A large number of metrics have been developed for this purpose, but there is currently no way of systematically testing the accuracy of individual metrics or knowing whether metrics provide consistent results. This lack of standardization can result in false-positive and false-negative findings if underpowered or inaccurate metrics are applied as tools for discovery. Here, we show that the choice of CUB metric impacts both the significance and measured effect sizes in numerous empirical datasets, raising questions about the generality of findings in published research. To bring about standardization, we developed a novel method to create synthetic protein-coding DNA sequences according to different models of codon usage. We use these benchmark sequences to identify the most accurate and robust metrics with regard to sequence length, GC content and amino acid heterogeneity. Finally, we show how our benchmark can aid the development of new metrics by providing feedback on its performance compared to the state of the art.

Diversity of Translation Initiation Mechanisms across Bacterial Species Is Driven by Environmental Conditions and Growth Demands.

The Shine-Dalgarno (SD) sequence motif is frequently found upstream of protein coding genes and is thought to be the dominant mechanism of translation initiation used by bacteria. Experimental studies have shown that the SD sequence facilitates start codon recognition and enhances translation initiation by directly interacting with the highly conserved anti-SD sequence on the 30S ribosomal subunit. However, the proportion of SD-led genes within a genome varies across species and the factors governing this variation in translation initiation mechanisms remain largely unknown. Here, we conduct a phylogenetically informed analysis and find that species capable of rapid growth contain a higher proportion of SD-led genes throughout their genomes. We show that SD sequence utilization covaries with a suite of genomic features that are important for efficient translation initiation and elongation. In addition to these endogenous genomic factors, we further show that exogenous environmental factors may influence the evolution of translation initiation mechanisms by finding that thermophilic species contain significantly more SD-led genes than mesophiles. Our results demonstrate that variation in translation initiation mechanisms across bacterial species is predictable and is a consequence of differential life-history strategies related to maximum growth rate and environmental-specific constraints.

Long-range correlations and fractal dynamics in C. elegans: Changes with aging and stress.

Reduced motor control is one of the most frequent features associated with aging and disease. Nonlinear and fractal analyses have proved to be useful in investigating human physiological alterations with age and disease. Similar findings have not been established for any of the model organisms typically studied by biologists, though. If the physiology of a simpler model organism displays the same characteristics, this fact would open a new research window on the control mechanisms that organisms use to regulate physiological processes during aging and stress. Here, we use a recently introduced animal-tracking technology to simultaneously follow tens of Caenorhabdits elegans for several hours and use tools from fractal physiology to quantitatively evaluate the effects of aging and temperature stress on nematode motility. Similar to human physiological signals, scaling analysis reveals long-range correlations in numerous motility variables, fractal properties in behavioral shifts, and fluctuation dynamics over a wide range of timescales. These properties change as a result of a superposition of age and stress-related adaptive mechanisms that regulate motility.

Leveraging genome-wide datasets to quantify the functional role of the anti-Shine-Dalgarno sequence in regulating translation efficiency.

Studies dating back to the 1970s established that sequence complementarity between the anti-Shine-Dalgarno (aSD) sequence on prokaryotic ribosomes and the 5' untranslated region of mRNAs helps to facilitate translation initiation. The optimal location of aSD sequence binding relative to the start codon, the full extents of the aSD sequence and the functional form of the relationship between aSD sequence complementarity and translation efficiency have not been fully resolved. Here, we investigate these relationships by leveraging the sequence diversity of endogenous genes and recently available genome-wide estimates of translation efficiency. We show that-after accounting for predicted mRNA structure-aSD sequence complementarity increases the translation of endogenous mRNAs by roughly 50%. Further, we observe that this relationship is nonlinear, with translation efficiency maximized for mRNAs with intermediate levels of aSD sequence complementarity. The mechanistic insights that we observe are highly robust: we find nearly identical results in multiple datasets spanning three distantly related bacteria. Further, we verify our main conclusions by re-analysing a controlled experimental dataset.

NullSeq: A Tool for Generating Random Coding Sequences with Desired Amino Acid and GC Contents.

The existence of over- and under-represented sequence motifs in genomes provides evidence of selective evolutionary pressures on biological mechanisms such as transcription, translation, ligand-substrate binding, and host immunity. In order to accurately identify motifs and other genome-scale patterns of interest, it is essential to be able to generate accurate null models that are appropriate for the sequences under study. While many tools have been developed to create random nucleotide sequences, protein coding sequences are subject to a unique set of constraints that complicates the process of generating appropriate null models. There are currently no tools available that allow users to create random coding sequences with specified amino acid composition and GC content for the purpose of hypothesis testing. Using the principle of maximum entropy, we developed a method that generates unbiased random sequences with pre-specified amino acid and GC content, which we have developed into a python package. Our method is the simplest way to obtain maximally unbiased random sequences that are subject to GC usage and primary amino acid sequence constraints. Furthermore, this approach can easily be expanded to create unbiased random sequences that incorporate more complicated constraints such as individual nucleotide usage or even di-nucleotide frequencies. The ability to generate correctly specified null models will allow researchers to accurately identify sequence motifs which will lead to a better understanding of biological processes as well as more effective engineering of biological systems.

A network approach to discerning the identities of C. elegans in a free moving population.

The study of C. elegans has led to ground-breaking discoveries in gene-function, neuronal circuits, and physiological responses. Subtle behavioral phenotypes, however, are often difficult to measure reproducibly. We have developed an experimental and computational infrastructure to simultaneously record and analyze the physical characteristics, movement, and social behaviors of dozens of interacting free-moving nematodes. Our algorithm implements a directed acyclic network that reconstructs the complex behavioral trajectories generated by individual C. elegans in a free moving population by chaining hundreds to thousands of short tracks into long contiguous trails. This technique allows for the high-throughput quantification of behavioral characteristics that require long-term observation of individual animals. The graphical interface we developed will enable researchers to uncover, in a reproducible manner, subtle time-dependent behavioral phenotypes that will allow dissection of the molecular mechanisms that give rise to organism-level behavior.

Depletion of Shine-Dalgarno Sequences Within Bacterial Coding Regions Is Expression Dependent.

Efficient and accurate protein synthesis is crucial for organismal survival in competitive environments. Translation efficiency (the number of proteins translated from a single mRNA in a given time period) is the combined result of differential translation initiation, elongation, and termination rates. Previous research identified the Shine-Dalgarno (SD) sequence as a modulator of translation initiation in bacterial genes, while codon usage biases are frequently implicated as a primary determinant of elongation rate variation. Recent studies have suggested that SD sequences within coding sequences may negatively affect translation elongation speed, but this claim remains controversial. Here, we present a metric to quantify the prevalence of SD sequences in coding regions. We analyze hundreds of bacterial genomes and find that the coding sequences of highly expressed genes systematically contain fewer SD sequences than expected, yielding a robust correlation between the normalized occurrence of SD sites and protein abundances across a range of bacterial taxa. We further show that depletion of SD sequences within ribosomal protein genes is correlated with organismal growth rates, supporting the hypothesis of strong selection against the presence of these sequences in coding regions and suggesting their association with translation efficiency in bacteria.

Dynamics and heterogeneity of a fate determinant during transition towards cell differentiation.

Yan is an ETS-domain transcription factor responsible for maintaining Drosophila eye cells in a multipotent state. Yan is at the core of a regulatory network that determines the time and place in which cells transit from multipotency to one of several differentiated lineages. Using a fluorescent reporter for Yan expression, we observed a biphasic distribution of Yan in multipotent cells, with a rapid inductive phase and slow decay phase. Transitions to various differentiated states occurred over the course of this dynamic process, suggesting that Yan expression level does not strongly determine cell potential. Consistent with this conclusion, perturbing Yan expression by varying gene dosage had no effect on cell fate transitions. However, we observed that as cells transited to differentiation, Yan expression became highly heterogeneous and this heterogeneity was transient. Signals received via the EGF Receptor were necessary for the transience in Yan noise since genetic loss caused sustained noise. Since these signals are essential for eye cells to differentiate, we suggest that dynamic heterogeneity of Yan is a necessary element of the transition process, and cell states are stabilized through noise reduction.

Impact of heterogeneity and socioeconomic factors on individual behavior in decentralized sharing ecosystems.

Tens of millions of individuals around the world use decentralized content distribution systems, a fact of growing social, economic, and technological importance. These sharing systems are poorly understood because, unlike in other technosocial systems, it is difficult to gather large-scale data about user behavior. Here, we investigate user activity patterns and the socioeconomic factors that could explain the behavior. Our analysis reveals that (i) the ecosystem is heterogeneous at several levels: content types are heterogeneous, users specialize in a few content types, and countries are heterogeneous in user profiles; and (ii) there is a strong correlation between socioeconomic indicators of a country and users behavior. Our findings open a research area on the dynamics of decentralized sharing ecosystems and the socioeconomic factors affecting them, and may have implications for the design of algorithms and for policymaking.

Changes in task-related functional connectivity across multiple spatial scales are related to reading performance.

Reading requires the interaction of a distributed set of cortical areas whose distinct patterns give rise to a wide range of individual skill. However, the nature of these neural interactions and their relation to reading performance are still poorly understood. Functional connectivity analyses of fMRI data can be used to characterize the nature of interactivity of distributed brain networks, yet most previous studies have focused on connectivity during task-free (i.e., "resting state") conditions. Here, we report new methods for assessing task-related functional connectivity using data-driven graph theoretical methods and describe how large-scale patterns of connectivity relate to individual variability in reading performance among children. We found that connectivity patterns of subjects performing a reading task could be decomposed hierarchically into multiple sub-networks, and we observed stronger long-range interaction between sub-networks in subjects with higher task accuracy. Additionally, we found a network of hub regions known to be critical to reading that displays increased short-range synchronization in higher accuracy subjects. These individual differences in task-related functional connectivity reveal that increased interaction between distant regions, coupled with selective local integration within key regions, is associated with better reading performance. Importantly, we show that task-related neuroimaging data contains far more information than usually extracted via standard univariate analyses--information that can meaningfully relate neural connectivity patterns to cognition and task.

Move-by-move dynamics of the advantage in chess matches reveals population-level learning of the game.

The complexity of chess matches has attracted broad interest since its invention. This complexity and the availability of large number of recorded matches make chess an ideal model systems for the study of population-level learning of a complex system. We systematically investigate the move-by-move dynamics of the white player's advantage from over seventy thousand high level chess matches spanning over 150 years. We find that the average advantage of the white player is positive and that it has been increasing over time. Currently, the average advantage of the white player is 0.17 pawns but it is exponentially approaching a value of 0.23 pawns with a characteristic time scale of 67 years. We also study the diffusion of the move dependence of the white player's advantage and find that it is non-Gaussian, has long-ranged anti-correlations and that after an initial period with no diffusion it becomes super-diffusive. We find that the duration of the non-diffusive period, corresponding to the opening stage of a match, is increasing in length and exponentially approaching a value of 15.6 moves with a characteristic time scale of 130 years. We interpret these two trends as a resulting from learning of the features of the game. Additionally, we find that the exponent [Formula: see text] characterizing the super-diffusive regime is increasing toward a value of 1.9, close to the ballistic regime. We suggest that this trend is due to the increased broadening of the range of abilities of chess players participating in major tournaments.